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Abstract 

We reexamine the conventional physical description of the neutrino 
evolution inside the Sun. We point out that the traditional resonance 
condition has physical meaning only in the limit of small values of 
the neutrino mixing angle, 8 -C 1. For large values of 9, the reso- 
nance condition specifies neither the point of the maximal violation 
of adiabaticity in the nonadiabatic case, nor the point where the fla- 
vor conversion occurs at the maximal rate in the adiabatic case. The 
corresponding correct conditions, valid for all values of 9 including 
9 > 7r/4, are presented. An adiabaticity condition valid for all values 
of 9 is also described. The results of accurate numerical computa- 
tions of the level jumping probability in the Sun are presented. These 
calculations cover a wide range of Am?, from the vacuum oscillation 
region to the region where the standard exponential approximation 
is good. A convenient empirical parametrization of these results in 
terms of elementary functions is given. The matter effects in the so- 
called "quasi- vacuum oscillation regime" are discussed. Finally, it is 
shown how the known analytical results for the exponential, 1/x, and 
linear matter distributions can be simply obtained from the formula 
for the hyperbolic tangent profile. An explicit formula for the jumping 
probability for the distribution N e oc (coth(x/Z) ± 1) is obtained. 



*New address from Sept. 1, 2000. 



1 Introduction 



The solar neutrino problem (SNP) is a discrepancy between the measured 
values of the solar neutrino flux at different energies 0J & H fl> HI and the 
corresponding predictions of the Standard Solar Model (SSM) |J. Not only 
is the observed flux suppressed, compared to the SSM predictions, but, if the 
data from the Homestake experiment are correct, the degree of suppression 
varies with energy. The leading explanation of this phenomenon is that 
neutrinos have small masses and the mass and flavor bases in the lepton 
sector are not aligned, just like in the quark sector. The resulting neutrino 
oscillations convert some of the solar electron neutrinos into another neutrino 
species. 

Neutrino oscillation solutions to the SNP have traditionally been divided 
into the so-called Mikheyev-Smirnov-Wolfenstein (MSW) solutions |7|, || |9| 
and the vacuum oscillation (VO) solutions, according to the physical mecha- 
nism responsible for the neutrino flavor conversion in each case. In the MSW 
case the conversion is caused by neutrino interactions with the solar (and 
Earth's) matter, while in the VO case it is due to long-wavelength neutrino 
oscillations in vacuum between the Sun and the Earth. Over time, it has 
become a tradition to treat the two cases completely separately, showing re- 
sults in separate plots (see, for example, |K], pT| ) and using different input 
formulas and different codes. 

Justifying such a complete separation, however, requires a careful analysis 
of the magnitude of the solar matter effects and the degree of decoherence of 
vacuum oscillations. The separation assumption has been recently reexam- 
ined by the author |12| and it has been found that the solar matter effects are 
nonnegligible for the vacuum oscillation solutions with Am 2 > 5x 10~ 10 eV 2 . 



This conclusion has been subsequently verified by other authors [13, 14, 15 



and the term quasi-vacuum oscillations (QVO) has been coined to refer to 
the region where both effects influence the neutrino survival probability (13 



It must be mentioned that the experimental situation has changed since 
the QVO solutions were first introduced. At the time, the most preferred 
part of the VO solutions was in the region Am 2 < 10~ 10 eV 2 . The latest 
Super-Kamiokande spectrum data, however, disfavors a large fraction of the 
vacuum oscillation region, roughly 2 x 10~ n eV 2 < Am 2 < 4 x 10~ 10 eV 2 
|T6| , pm . At the same time, the solutions with Am 2 > 4 x 10" 10 eV 2 , i. e. 
the QVO solutions, remain allowed. 

Prior to |T3, the VO solutions had always been studied in the range of the 
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neutrino mixing angle < 9 < 71/ 4 for a fixed sign of Am 2 . When matter 
effects are included, however, this only covers a half of the full parameter 
space. To cover the full space, one can either (i) keep 9 in the range < 9 < 
7r/4 and consider both signs of Am 2 , or (ii) fix the sign of Am 2 and vary 9 
from to 7r/2. We advocate the second option as a better physical choice, 
because it makes manifest the continuity of physics around the maximal 
mixing ]17|, fjj| . 

The parametrization < 9 < ir/2 requires one to reexamine the choice 
of a variable for plots, because the traditional variable sin 2 29 is not suitable 
for this purpose JD]. While either 9 or sin 2 9 would be adequate for plotting 
only the QVO region, neither choice allows one to take a global view of the 
neutrino parameter space and show all solutions, including the SMA and 
(quasi) vacuum oscillation solutions, on the same plot |i~8| . A particularly 



convenient choice turns out to be tan 9 on the logarithmic scale, first used 
in |19] to describe 3-family MSW oscillations. In addition to covering the 
range < 9 < 71/ 2, it also does not introduce any unphysical singularity 
around 9 = 7r/4 (unlike the traditional sin 2 29, see [|17|) and makes it easy to 
see where in the vacuum oscillation region the evolution in the Sun becomes 
completely nonadiabatic (points 9 and tt/2 — 9 become equivalent, so that 
solutions become symmetric with respect to the 9 = tt/A line). 

In the first part of this paper we address several conceptual questions 
that arise in the analysis of the solar matter effects and become particularly 
apparent for the values of the mixing angle 9 > 71 /A. To introduce these 
questions, it is useful to first summarize the basic mechanism of the neu- 
trino evolution in the Sun. Inside the Sun, because of the changing electron 
density, the eigenstates of the instantaneous Hamiltonian change along the 
neutrino trajectory. For Am 2 /E u <C 10~ 5 eV 2 /MeV the neutrino is produced 
almost completely in the heavy eigenstate. If the parameters Am 2 /E„ and 
9 are such that the neutrino remains in the heavy eigenstate as it travels to 
the solar surface (adiabatic evolution), there are no subsequent vacuum oscil- 
lations. To oscillate in vacuum, as a necessary condition, the neutrino must 
at some radius in the Sun "jump" into the superposition of the heavy (i/ 2 ) 
and light (vi) mass eigenstates. Conventional VO regime is realized when 
this "jumping" is extremely nonadiabatic (preserving flavor), in which case 
the neutrino exits the Sun as cos9\v\) + e 1 ^ sin 6*1^2). In the QVO regime 
the neutrino still partially jumps in the v\ eigenstate, but with a smaller 
amplitude. 

The obvious questions one would like to answer are: 
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• What physical criteria determine whether the neutrino evolution is adi- 
abatic or not? 

• At what radius in the Sun does the nonadiabatic "jumping" between 
the eigenstates of the instantaneous Hamiltonian take place? 

The traditional wisdom is that one should analyze the density profile around 
the so-called resonance point, i.e., the point where the difference of the eigen- 
values of the instantaneous Hamiltonian is minimal and the local value of the 



mixing angle is 6 m = ir/A (see, e.g., |2(| [21], |3|, |4[ ^5] |26|). This, how- 



ever, clearly needs to be modified for large mixing angles. In particular, for 
6 > 7r/4 the resonance, defined in this way, simply does not exist. We will 
show how this contradiction is resolved in Section 12. 2L In Section 12.31 we for- 



mulate the adiabaticity condition that, unlike the standard result, remains 
valid for 9 > n/4. 

In Section |^ we present the results of numerical calculations of the jump- 
ing probability P c for the neutrino propagating in the realistic solar profile. 
The calculations are carried out for a wide range of Am 2 and tan 2 6, from 
the VO region to the region where the exponential density approximation is 
valid. We show how the adiabaticity prescription of Section |2]3| applies to 
this case. We also give a simple empirical prescription on how to compute P c 
in this range of the parameters in terms of only elementary functions. Such 
an empirical parametrization of the numerical results can be helpful if one 
would like to be able to quickly estimate the value of P c anywhere in the 
range in question without having to solve the differential equation each time. 



In Section |3.2| , we discuss what happens in the transitional region between 
the adiabatic and nonadiabatic regimes (QVO). In particular, we determine 
what part of the solar electron density profile is primarily responsible for the 
matter effects in this region. 

Finally, in Section |4] we comment on the four known exact analytical 
solutions for the neutrino jumping probability P c . Such solutions have been 
found for the linear, exponential, 1/r, and the hyperbolic tangent matter 
density profiles. A natural question to ask is whether these profiles have 
something in common that makes finding exact solutions possible. Using 
the formulation of the evolution equations introduced in Section we show 
that all four results are not independent and that, given the formula for 
the hyperbolic tangent profile, one can very simply obtain the other three 
solutions. As an added benefit, we obtain an exact expression for the density 
distribution N e oc (coth(x/l) ± 1). 
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2 Physics of the nonadiabatic neutrino evo- 
lution 

2.1 Review of the oscillation formalism 

For completeness, we begin by summarizing the well-known basic formalism 
for neutrino oscillations in matter. In the simplest case, when the mixing is 
between v e and another active neutrino species, the evolution of the neutrino 
state is determined by four parameters: the mass-squared splitting between 
the neutrino mass eigenstates Am 2 = m\ — m\ , the neutrino mixing angle 9, 
the neutrino energy E u , and the electron number density N e of the medium. 
One has to solve the Schrodinger equation id<fi/dt = H<p, where = (<p e , (j)^) T 
is the state vector made up of the electron neutrino and the muon neutrino.[| 
The Hamiltonian H is given by [0 



where A = Am 2 / (4^) and A = \ / 2G F N e /2. The constant in the Hamilto- 
nian is irrelevant for the study of oscillations and will be omitted from now 
on. The time variable t may be replaced by the distance traveled x, since 
the solar neutrinos are ultrarelativistic. 

For a constant electron number density iV e the Hamiltonian can be triv- 
ially diagonalized, H' = VHV* = diag(— A m , +A m ), where 

A m = \J (A - A cos 29) 2 + (A sin 29) 2 = VA 2 -2AA cos 29 + A 2 . (2.2) 
In terms of A m , the Hamiltonian ( |2.1| ) can be rewritten as 



where 9 m is the mixing angle in matter. The rotation matrix V is given by 



*In reality, </> M here denotes a linear combination of cj)^ and <f> T in which <j) e oscillates. 




(2.1) 




(2.3) 




(2.4) 
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The parameters A m and 6 m are related to the original parameters in the 
Hamiltonian ( |2.1|) as follows 

A m sin20 m = A sin 26, (2.5) 
A m cos2# m = A cos 26 -A, (2.6) 

A sin 26 , , 

tan2# m = -—— . 2.7 

A cos 26 — A 

We will always label the light mass eigenstate by v\ and the heavy one 
by u 2 . Since one can redefine the phases of v e ^ and v\^-, it is easy to see 
that in this convention the physical range of the mixing angle is0<#<7r/2 



18, 17 



As long as N e (x) is constant, the time evolution of the mass eigenstates 
is particularly simple. Each of the two mass eigenstates evolves only by 
a phase: \vi(t)) = |^i(0)) exp(iA m t), l^it)) = [^(O)) exp(— iA m t). If at 
time t = the neutrino state is a linear combination a|^i(0)) + 6|z/ 2 (0)), the 
absolute values of the coefficients a and b do not change with time, i.e., the 
probability for the neutrino to "jump" from one Hamiltonian eigenstate to 
another, P c = \a(t = +oo)| 2 — \a(t = 0)| 2 , is trivially zero. 

Consider next the case of a varying electron density. In this case, in 



general, one can no longer diagonalize the Hamiltonian in Eq. (2.1). How- 
ever, one can still speak of the eigenstates of the instantaneous Hamiltonian 
(henceforth, "the matter mass eigenstates"), and define the jumping proba- 
bility between those states. It turns out that if the electron density changes 
sufficiently slowly along the neutrino trajectory (to be quantified later), the 
jumping probability vanishes, just like in the constant density case. This 
is known as the adiabatic evolution. At the same time, when the density 
changes abruptly, the jumping probability is clearly nonzero. In particular, 
if the neutrino crosses a step-function density discontinuity, the flavor state 
does not have any time to evolve, while the mass basis in matter instanta- 
neously rotates. It is easy to see that in this situation, known as the extreme 
nonadiabatic evolution, the jumping probability is given by 

P C NA = sin 2 (fl bcforc - # after ). (2.8) 

In general, for a monotonically varying density P c lies between and P C NA . 

In this paper we are concerned with the evolution of solar neutrinos. The 
electron number density inside the Sun, N e , falls off as a function of the 
distance from the center r as shown in Fig. 0. In the range 0.15_R Q < 
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Figure 1: The electron number density profile of the Sun according to the 
BP2000 standard solar model. 



r < 0.65-Rq the profile can be approximated very well by an exponential 
N e (r) cx exp(— r/r ), with r = -R Q /10.54 = 6.60 x 10 4 km (shown by a 
straight line in the Figure). However, in the convective zone of the Sun and 
also in the core where the neutrinos are produced the profile deviates rather 
significantly from exponential. 

In order to study the jumping probability between the matter mass eigen- 
states it is convenient to change to the basis these states define. Substituting 
in the Schrodinger equation (ft = V^ip, ip = (/0i, 1P2), we get 

i-^-(Vty) = iZVty, 
ax 

dip , dV^ , s 

i-r- = VHV^ip -iV—ip. 2.9 
dx dx 
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Since VHV^ = diag(— A m , +A m ) and from Eq. 

dV* ( 1 \ da 



in 



dx \ — 1 J dx ' 
we obtain the desired evolution equation in the basis of the matter mass 



eigenstates |20 



— ( ^ \ = ( -d^ra/dx \ ( Ipi 

dx\ ^2 I \ dOm/dx -iA m I \ ip 2 



(2.10) 



The steps outlined so far are standard in the treatment of the MSW 
effect. We will next make an extra step that will prove very helpful for 
the subsequent analysis, particularly in Section f|. Namely, we will choose 
9 m instead of independent variable. So long as the density varies 

monotonically, such a change is one-to-one. Eq. ( |2.10|) becomes 



d ( ipi\( iA m /e m -l \( 

d9 m \^2 { 1 -iA m /9 m [ip 2 



(2.11) 



Here A m and 9 m = d9 m /dx can both be expressed in terms of 9 m using the 
following relationships 

• sin 2 2g m dA 

m ~ 2 A sin 20 da;' [ ' 

\n 2A 2 sin 2 2fl 1 
9 m sm 3 29 m dAjdx 

A = Asin ( 2f) -- 29 », (2.14) 
sin2fl m ' v ; 



which follow directly from Eqs. (|2.5| - |2.7| ). For instance, for the infinitely ex- 
tending exponential profile A(x) = A exp(— x/r ) the derivative is dA(x) /dx = 
— A(x)/r and so 

Am = 2Ar sin 2 2fl 

9 m sm 2 29 m sm(29 m -29)' 1 ' ' 

The angle 9 m varies from its value at the production point 9 & to its 
vacuum value 9. For the infinite exponential profile we have 9 Q — > ir/2. 
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Notice that the quantity A m /9 m in Eq. (|2.15|) is singular when 6 m approaches 
either of its limiting values, as should be expected. 

The shape of the function |A m /0 m | for 6 = tt/4 and A = 10" 9 eV 2 /MeV 
for the idealized exponential profile is shown in Fig. |2| by the solid curve. The 
value of r was chosen to be i? o /10.54, the slope of the best fit line in Fig. [I]. 
The dashed curve in the Figure shows the same quantity for the realistic 
BP2000 solar profile for the same values of 9 and A. It is important to keep 
in mind that the two curves change qualitatively differently as one changes 
A. While the exponential curve just scales by an overall factor, the BP2000 
curve also changes its shape, approaching the shape of the (rescaled) solid 
curve for sufficiently large values of A. 
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2.2 Modification of the notion of resonance for large 0. 

We are ready to address the questions posed in the introduction. A con- 
venient starting point is the evolution equation in the form of Eq. (p.ll|) . 
Neutrino is produced in the v e state, which corresponds to (cos 9 Q , sin 9 Q ) in 
the matter mass basis. To simplify the presentation, we shall consider the 
values of A 10~ 5 eV 2 /MeV, so that at the production point A <C A Q , and 
hence (cos 9 Q , sin 9 & ) — > (0,1). 

The evolution equation can be easily solved in the two limiting cases. If 
the off-diagonal elements can be neglected, the evolution is adiabatic, i.e., 
\ip2{@m)\ 2 remains constant. If, on the other hand, for almost the entire 
interval between 9 & and 9 the diagonal terms can be neglected, the solution 
is 

ip{9) = exp 

l J o@ \ 1 u / J 

(2.16) 

This limit corresponds to the extreme nonadiabatic case. The corresponding 
jumping probability equals P c = sin 2 (^ — 9) = cos 2 9, in agreement with 
Eq. (p|). 

Returning for a moment to the physical x-space, we note that no jumping 
between the mass eigenstates occurs either in the solar core [[jj or in vacuum. 
The nonadiabatic evolution takes place in a localized region, with a center 
at the point of "the maximal violation of adiabaticity" . Our goal next is to 
establish the location of this point. 

Conventional wisdom says that the adiabaticity condition is violated max- 
imally at the resonance point 

A = A cos 20, (2.17) 

where the separation between the energy levels is minimal and 9 m = 7r/4. 
This assertion can be found in the early papers PTl |27j[], as well as in 

numerous subsequent reviews on the subject, e.g. [|4], [25], |26|. However, in 
all these papers it is assumed either explicitly or implicitly that the vacuum 




'''A notable exception is Ref. [2£]. We do not agree, however, with the adiabaticity 



criterion proposed there (see Section 2.3) 
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mixing angle 9 is small. It is easy to see that for a large value of the mixing 
angle the use of the condition in Eq. (|2.17| ) leads to a contradiction. 

For small 9, Eq. Q2.17| ) is satisfied in a layer in the Sun where the density 



is A(x) ~ A. As the value of 9 increases, Eq. (|2.17|) predicts that the 
resonance occurs at lower and lower electron density until, as 9 approaches 
7r/4, it moves off to infinity. It is not obvious how to interpret the last result, 
as it is physically clear that no level jumping can occur at infinity where 
neutrinos undergo ordinary vacuum oscillations. The difficulty is even more 
obvious when 9 > 7r/4, in which case the resonance simply never occurs. At 
the same time, as already mentioned, in the extreme nonadiabatic regime the 
level jumping probability is nonzero for any value of 9 and varies smoothly 
around 9 = tt/4, P c = cos 2 9. 

The resolution to this apparent paradox is very simple. As Eq. (|2.11|) 
suggests, adiabaticity is maximally violated at the minimum of \A m /9 m \ 
( [p0| , p7|). It is easy to see, however, that the minimum of \A m /9 m \ in general 
does not reduce to the condition of Eq. ( p. 17] ) . This can be explicitly seen on 
the example of the infinite exponential density distribution. Differentiating 
Eq. (|2.15|) , one finds that the minimum occurs when 

cot(20 m - 29) + 2 cot(2# m ) = 0, (2.18) 

or 



Unlike Eq. ( gT7| ), Eq. fl2jj ) states there is a nonadiabatic part of the 



neutrino trajectory for all physical values of 9, including 9 > ir/4. While both 
equations for small 9 predict that jumping between the local mass eigenstates 
occurs around A = A, Eq. ( [2.19 ) states that for maximal mixing it occurs 



around A = A/ v2, not at infinity, and for 9 close to 7r/2 it happens around 
A = A/2, all physically sensible results. 

The situation is illustrated in Fig. |3|, which shows the probability of 
finding the neutrino in the heavy mass state z/2 as a function of the distance 
x. The parameters of the exponential were taken from the fit line in Fig. p 
and Am 2 /E u = 10~ 9 eV 2 /MeV. Three large values of the mixing angle (9 = 
7r/6, 7r/4, and 7r/3) and one small value {9 = 7r/60) were chosen. The dashed 
lines and dots mark the points where adiabaticity is maximally violated, as 
predicted by Eq. ( |2.19| ). One can see that the partial jumping into the light 
mass eigenstate in all four cases indeed occurs around the marked points. 
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1^2 




Figure 3: Neutrino state evolution in the case of the infinite exponential 
density profile for Am 2 /E„ = 10~ 9 eV 2 /MeV. The plot shows the probability 
of finding the neutrino in the heavy matter mass eigenstate u 2 as a function 
of position x, for four different values of the vacuum mixing angle. The 
points of the maximal violation of adiabaticity, as predicted by Eq. (|2.19 ), 
are marked. 



It is instructive to analyze each of the factors A m and l/9 m separately. 
The first one indeed has a minimum at the traditional resonance point 9 m = 
7r/4 (corresponding to Eq. fl2.17|) ) or, if 9 > 7r/4, at the endpoint 9 m = 9. The 
second one, however, has a minimum at a point which is, in general, different 
from the resonance. This minimum exists for all values of 9, including 9 > 
ir/ 4. In the case of the exponential profile, it is located at 

9 m = vr/4 + 9/2, (2.20) 

or halfway between 9 and 9 Q = tt/2. The corresponding value of the density 
at that point is 

A = A. (2.21) 

This coincides with the resonance condition for ^ C 1 and that is why the 
standard resonance description works very well in this limit. In general, 
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however, the minimum of the product function lies somewhere between 7r/4 
and vr/4 + 9/2 (or 9 and vr/4 + 9/2, if 9 > vr/4)|. 

It is also worth mentioning that Eqs. ( p.20| ) and fl2.21| ) represent an 
important condition in the case of the adiabatic evolution. Namely, they 
specify a point where the rate of rotation of the mass basis with respect to 
the flavor basis is maximal, which in the adiabatic case can be interpreted 
as a point where the flavor composition of the neutrino state changes at the 
fastest rate. This shows that for large 9 Eq. ( |2.17|) not only does not describe 
the point of the maximal violation of adiabaticity in the nonadiabatic regime, 
but also does not specify the point where the flavor conversion occurs at the 
maximal rate in the adiabatic regime. 

To summarize, Eq. ( |2.17| ) can only be used in the small angle limit. Even 
in that case one should be careful applying it for certain purposes. Note, for 
example, that Eqs. ( |2.17| ) and ( |2.19| ) have different Taylor series expansion 
around 9 = 0, 



A 
A 



A(l 
A(l 



29 2 /3) 
29 2 ) 



for Eq. (2.19), 
for Eq. ([247]). 



Thus, even at small 9, Eq. (p. 17 ) fails to predict how the point of maximal 
nonadiabaticity shifts function of 9. 

The belief that jumping between the matter mass eigenstates occurs at 
the resonance for all values of 9 might have been one of the reasons for the 
tradition to treat separately the cases of 9 < 7r/4 and 9 > 7r/4, obscuring 
the fact that physics is completely continuous across 9 = 7r/4. Over the 
years, it has caused some unfortunate confusions, as exemplified by the flawed 
criticism of the results of [0 in 29fl . It was probably the principal reason 



why the correct expression for the electron neutrino survival probability in 
the 9 > 7r/4 part of the QVO region was not given until recently |17|, [L2| (c./. 
Eq. (6) in @). 



One important application of Eq. ( [2.19 ) is the determination of the phase 
of vacuum oscillations on the Earth |3{], fJT]. If jumping between matter mass 
eigenstates indeed occurred at the resonance ( P-ITQ , one would expect that 
in the canonical vacuum oscillation formula 



P = 1 - sin 2 29 sin 2 ( 1.27^-^ + 5 T 
\ E 



[2.22) 



Ht is even possible for certain density profiles and certain values of A and 9 to have 
more than one minima. 
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the residual phase <5 res would be minimized when the distance L was measured 



from the resonance in the Sun. Ref. [[H] indeed begins with this assumption, 
but after presenting the resulting formulas notes that the residual phase is 
much smaller if L is instead measured from the layer where A = A, not 
A = A cos 29. It is unfortunate that this important observation has not 
received proper attention and was not further developed in the subsequent 
literature. The preceding discussion shows that this result is not just a 
mathematical coincidence, but has a simple physical explanation. 

Finally, it is important to discuss at what density the adiabaticity is 
maximally violated in the case of the realistic solar profile. Qualitatively, it 
is not difficult to anticipate the changes to Eq. fl2.19|) in this case. 



For sufficiently large A, the adiabaticity is maximally violated in the 
radiative zone, entirely within the exponential part of the profile 
that Eq. ( ^.1^ ) directly applies^ 



so 



For small A, the nonadiabatic part of the trajectory lies close to the 
surface of the Sun where the profile falls off rather rapidly. While 
for small 9 the minimum of A m /9 m should still, of course, occur at 
A = A, for large angles it is shifted to values of A somewhat lower 
than those predicted by Eq. ( |2.19| ). The evolution in the latter case 



will be discussed in more detail in Section 3.2 



• For A in the intermediate range, the jumping occurs near the bottom 
of the convective zone where the density falls off somewhat slower than 
in the exponential part. There the value of the ratio A/ A at large 9 
should increase compared to Eq. ( |2.19| ). 

These qualitative expectations are supported by the results of numerical 
calculations presented in Fig. |], where the ratio Aj A at the point of minimal 
A m /(9 m is plotted as function of 9. The three curves shown correspond to the 
values of A in the three different regimes: A = 2 x 10~ 6 eV 2 /MeV (curve 1), 
A = 7 x 1(T 8 eV 2 /MeV (curve 2), and A = 1 x 1CT 8 eV 2 /MeV (curve 3). 
The numerical studies of the BP2000 solar profile are presented in Section |3|. 



§As we shall see in Section 
for 6 -C 1, in which case Eq. (2.19) reduces to A 



3, in this range of A the nonadiabatic evolution only occurs 
A. 
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o [ 

0.25 0.5 0.75 1 1.25 1.5 



Figure 4: The value of the ratio A(x) /A at the minimum of A m / 8 m computed 
for the realistic solar density profile (BP2000). The curves correspond to 
A = 2 x 10" 6 eV 2 /MeV (curve 1), A = 7 x 10" 8 eV 2 /MeV (curve 2), and 
A = 1 x lO" 8 eV 2 /MeV (curve 3). 



2.3 Adiabaticity condition for large 9. 

We now turn to formulating the adiabaticity condition that is valid for all, 
and not just small, values of the mixing angle 9. At first sight it appears 
that for the evolution to be adiabatic it is simply enough to require that the 
diagonal elements in Eq. fl2.ll ) be larger than the off-diagonal ones. This 



condition becomes most critical at the point where the adiabaticity is max- 
imally violated. Since traditionally this point has been identified with the 
resonance, a commonly cited condition is [p0[ 



|A m /# m |<9 m=7r / 4 > 1. (2.23) 

Since we have shown that the point of the maximal violation of adiabatic- 
ity in general does not coincide with the resonance, Eq. fl2.23p clearly needs 
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to be modified for large mixing angles. Superficially, the problem appears 
easy to fix: one should evaluate the left hand side at a value of 9 m close to 
the true point of the maximal violation of adiabaticity. In the case of the 
exponential profile, this point would be given by the solution of Eq. (|2.18|) . 
For the purpose of an estimate, we can approximate it by the point halfway 
between 9 and vr/2 (see Eq. (gggfl), 

> 1. (2.24) 

This condition, however, still turns out to be inadequate for large values 
of 9. To obtain the correct condition, one has to analyze the problem more 
carefully. 

The key is to express the information contained in the system of two 
evolution equations in a single equation. Such an equation can be easily 
written for the ratio s = ipi/ip2- From Eq. (|2.11|) it follows that s obeys the 
following first order equation 

(ls 2i^s-(s 2 + l). (2.25) 



d9 m 9 m 

It is easy to see that by neglecting appropriate terms on the right hand 
side one obtains both the adiabatic and nonadiabatic limits. The adiabatic 
limit corresponds to neglecting the terms in parentheses, while the extreme 
nonadiabatic limit is obtained if one neglects the first term on the right. 
In the second case the solution is s = cot(9 m ). Thus, the self consistent 
condition to have the extreme nonadiabatic solution is 2| A m /9 m s\ <C (s 2 + 1), 
with s = cot(6 l m ), 9 m = 7r/4 + 9/2. In the opposite limit, the evolution is 
adiabatic. Thus, we obtain the adiabaticity condition 

\A m /9 m \ dm ^ /4+e/ 2 > (tan(7r/4 + 9/2) + cot(7r/4 + 9/2))/2, (2.26) 

or, upon simplification, 

cos9\A m /9 m \o m ^ 7r / 4:+ e/2 > 1. (2.27) 

Notice, that all three conditions, Eqs. (|2T23|) , flggj ), and (|2~2p agree in 



the small angle limit. For large angles, however, especially for 9 > tt/4, they 
give very different predictions and only Eq. ( |2.27|) is the correct one. 

Let us demonstrate this on the example of the infinite exponential profile. 
In this case the exact analytical solution is known [32 , |33| (see also |34fl ). 



7Cos 2 6» _ ^ 

P c = . (2.28) 
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Figure 5: The boundary between the adiabatic and nonadiabatic regions in 
the case of the exponential profile, as predicted by Eq. (|2.23| ), Eq. ( |2.24| ) , 
and Eq. (|2.27|) . The dashed curves show the contours of constant P c . For 
large angles Eq. fl2.27j ) provides a better description of the boundary. 
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Eq. ( |2.28| ) has two regimes. For small mixing angles, the formula reduces to 
P c = exp(— 7 sin 2 9), so the evolution is nonadiabatic when 76 12 = 47ir A9 2 <C 
1. By contract, for 9 > 71/ 4 and 7 ^> 1 the numerator is always much 
smaller than the denominator. The transition between the adiabatic and 
nonadiabatic regimes now occurs when 7 ~ 1, with a weak dependence on 9. 

Let us see if Eq. ( |2.27| ) correctly captures this behavior. Using Eq. ( |2.15| ) 
we obtain 

8Ar sin 2 #>l. (2.29) 

Fig. |5] shows the contour of SArosin 2 ^ = 1 computed for r = -R Q /10.54 
(solid line). For comparison, the dash-dotted curves shows the corresponding 
prediction of Eqs. ( 2.23 ) and ( p. 24 ). The dashed curves are the contours of 



constant P c computed using Eq. ( |2.28| ). It is clear from the Figure that the 



description of Eq. (|2.27|) is correct not only for small 9, but also for 9 > 7r/4. 
The transition from adiabaticity to nonadiabaticity for 9 > tt/4 occurs for 
Am 2 /E u ~ 1CT 9 eV 2 /MeV, precisely where 7 ~ 1. 

The true usefulness of Eq. fl2.27|) is not so much in being able to explain 



the physics behind the known analytical solution as it is in being able to make 
predictions for a variety of new density profiles, provided those profiles are 
sufficiently smooth. Such an analysis, however, would be beyond the scope of 
this paper. The question we do want to address is what Eq. ( |2.27| ) predicts for 
the realistic solar density profile. It turns out that for large mixing angles and 
the values of A that yield nonadiabatic evolution the nonadiabatic jumping 
takes place mostly in the convective zone, where the density profile deviates 
from the exponential, and so one can no longer rely on ( |2.28| ) to get the 
shape of the boundary between the adiabatic and nonadiabatic regions. Eq. 
( |2.27| ) nonetheless works quite well. The corresponding numerical results are 
presented next. 



3 Calculations for the realistic solar profile 

3.1 Numerical calculations of the jumping probability 

In this Section we present the results of numerical computation of the jump- 
ing probability P c with the realistic (BP2000) solar density profile. In prac- 
tice, ultimately, one would like to compute expected event rates at various 
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experiments, and for that one needs to know the solar neutrino survival prob- 
ability P{v e — > u e ). It is, however, quite straightforward to show that, for a 
given neutrino energy, the survival probability is given by 

P(v e -> u e ) = Pi cos 2 9 + (1 - Pi) sin 2 9 



/ Am 2 \ 

P c (l-P c )cos2fl sin2#cosl2.54-^L + <5j , (3.1) 



where Pi = P c sin 9 & + (1 — P c ) cos 9® [0. Thus, the problem of finding 
P{y e — > u e ) reduces to finding the jumping probability P c . 

The quantity P c can be found analytically in several limiting regimes. For 
Am 2 /E u > 10~ 7 eV 2 /MeV the condition of Eq. ( |2.19| ) is satisfied well inside 
the Sun where the profile is exponential with tq = P Q /10.54. Hence, in this 
case the jumping probability should be adequately described by Eq. (|2.28| ). 



The second regime is 9 <C 1, in which case the standard resonance de- 
scription applies and, importantly, the resonance is very narrow. As a result, 
the jumping probability can be adequately described by the analytical for- 
mula for a linear density profile, P c = exp(— tiA 2 sin 2 29\dA(x) / dx\^ =A ) (Eq. 



Q4.1J) of Sect. |]). The contours of constant P c in this regime are expected 
to follow the behavior of the changing slope of the BP2000 density profile 
shown in Fig. [I]. 

The third regime is the region of small A and large 9. As already dis- 
cussed, for A — > the evolution becomes extremely nonadiabatic, i. e., 
P c — > cos 2 9. For Am 2 /E u ~ 10~ 9 eV 2 /MeV we have the so-called quasivac- 
uum oscillation region. It was shown in |12|] that for Am 2 / E v £ 5 x 10~ 9 
eV 2 /MeV Eq. (|2T28| ) provides a good fit to numerical calculations, provided 



one takes ro = Plq/18.4. This region will be discussed in more detail in 



Section 3.2 



We next present numerical results for P c that cover the entire range be- 
tween these regimes. We numerically solve Eq. (|2.25|) on a grid of points in 
the range 1(T 3 < tan 2 # < 10, l(T n eV 2 /MeV < Am 2 / E v < 2xl0" 7 eV 2 /MeV. 
Notice that, while the same result would be obtained by solving the system of 
equations in Eq. ( |2.1(J| ), Eq. fl2.25p requires significantly less computer time. 



Indeed, Eq. ( |2.10| ) contains four real functions (real and imaginary parts of 
ipi and ^2), but only two of them are independent, since |-?/>i| 2 + |^ 2 | 2 = 1 
and the overall phase has no physical meaning. 

The resulting contours of constant P c are shown in Fig. ^ by solid curves. 
All the features anticipated in the discussion above are clearly present. 
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We can now test the validity of the adiabaticity criterion introduced in 
the previous Section (Eq. ( |2.27| )). The shaded region in Fig. [I] corresponds 
to cos6 l |A m /^ m |e m=7r /4 + 0/2 > 1. As one can see by comparing its shape with 
that of the contours of constant P c , the criterion in question indeed describes 
the adiabatic region quite well, even for 9 > 7r/4. 

For practical calculations it is convenient to have a simple expression 
involving only elementary functions that provides a satisfactory means of 
estimating P c without having to solve the differential equation each time. 
We can construct such a purely empirical fit function by taking Eq. ([2.28Q as 
a starting point. Since for both large and small values of Am 2 / E v one can 
use Eq. (|2.28|) with the appropriate values of r , we can modify Eq. (|2.28|) by 



making r a function of Am 2 /E u , smoothly interpolating between P /1O.54 
and #0/18.4, 



r Q {Am 2 /E v 



Re 
10.54 



0.75 ( - arctan[-10(a + 8)1 + - ) + 1 

V 7T 2 ' 



-1 



(3.2) 



where a = log 10 [(Am 2 /Pj,)/(eV 2 /MeV)]. This provides an adequate fit for 
both Am 2 /E u < 5 x 10~ 9 eV 2 /MeV and Am 2 / E u > 10~ 7 eV 2 /MeV. By 
making additional modifications to the function it is possible to obtain a 
reasonable fit also to the contours in the intermediate region, 



rl\Am 2 /E v ) 




1 i 
0.75 ( - arctan[-10(a + 7.95) + 0.96] + - 

7T 2 



0.5 2 i 



n -1 



(3.3) 



where b = log 10 [tan 2 9}. It is important to emphasize that this equation 
should only be viewed as a purely empirical fit to the numerical results, 
intended to facilitate practical computations of the solar neutrino survival 
probability. 

The contours of constant P c fit computed with Eq. (|3.3| ) are shown by 
dashed lines in Fig. |6|. The discrepancy between the numerical results and 
P c fit is |P c fit -P c BP2000 | < 0.022. 

3.2 Matter effects in the QVO region. 

In this subsection we discuss the matter effects in the QVO regime. This 
region, characterized by Am 2 /E u ~ 10~ 9 eV 2 /MeV and large values of the 
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Figure 6: Contours of constant level jumping probability P c . The solid line 
represent the result of numerical calculations using the BP2000 solar profile. 
The dashed lines illustrate a possible empirical fit using only elementary 
functions (Eq. (gjg)). 
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Figure 7: The comparison of the neutrino state evolution in the infinite 
exponential profile (dashed) and the realistic solar profile according to the 
BP2000 solar model (solid), plotted as a function of the mixing angle in 
matter 9 m . The bottom figure shows the deviation of the evolution in both 
cases from extreme nonadiabatic. 
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Figure 8: Same as Fig. [7|, but as a function of the distance from the center 
of the Sun x. 
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mixing angle 9 deserves a separate treatment. While lying below the adia- 
batic region as determined by the criterion of Eq. (|2.27| ) (see Fig. ||), it is 
nonetheless characterized by a significant deviation from extreme nonadia- 
baticity. In fact, as discussed in the last section, for Am 2 /E u < 5 x 10~ 9 
eV 2 /MeV the jumping probability in this region can be described quite well 
by Eq. ( |2.28| ) with r = -R /18.4. Thus, the neutrino evolution in this region 
is partially adiabatic, and our task next is to sketch a simple physical picture 
of this phenomenon. 

First, notice that the condition A ~ A in the QVO region is satisfied very 
close to the Sun's edge (at x = 6.6 x 10 5 km for Am 2 / E u = 10~ 9 eV 2 /MeV), 
where the profile falls off very rapidly. This of course does not mean that the 
evolution is extremely nonadiabatic. As discussed previously, the slope of the 
profile at the point A = A can only be used to find the jumping probability 
in the case of small mixing angles. For large mixing angles, the nonadiabatic 
segment of the neutrino trajectory is large, and roughly the first half of this 
segment lies within x < O.92i? = 6.35 x 10 5 km, where the density is equal 
to or greater than that of the fitted exponential profile. 

This suggests that the evolution of the solar neutrino in the QVO regime 
can be modeled as the evolution in the exponential profile truncated near the 
point of the maximal violation of adiabaticity. Let us investigate the basic 
features of this model. 

First, consider a partially adiabatic evolution (cos 9\ A m / 9 m \ 9m ^ n / 4:+e /2 ~ 
1) in case of the infinite exponential profile. As the neutrino traverses the 
Sun, its state vector, although not completely "attached" to the heavy mass 
eigenstate, rotates away from the pure v e flavor state. This rotation occurs 
both before and after the point of the maximal violation of adiabaticity. 
Furthermore, an important role is played by the parts of the trajectory where 
9 m is close to either it/ 2 or 9, where the quantity A m /9 m is large (see Fig. fj). 

Now truncate the profile close to the point of the maximal violation of 
adiabaticity. Then the amount of the rotation of the state with respect 
to the flavor basis decreases roughly by a factor of two. Actually, because 
the function in Eq. Q2.15D has a double pole at 9 m = 7r/2, and a single pole 
around 9 m = 9 the factor can be expected to be somewhat less than two. The 
numerical calculations of the last section indeed show that for the BP2000 
profile the parameter r$ changes from i? Q /10.54 to i? Q /18.4. 

To illustrate this more quantitatively, let us see how \ip2{9 m )\ deviates 
from its extreme nonadiabatic behavior \ip2{9 m )\ = sin9 m in the QVO regime. 
Fig. 0(A) shows the probability of finding the neutrino in the heavy mass 
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eigenstate, for Am 2 /E u = 10~ 9 eV 2 /MeV, 9 = 7r/4, as a function of the 
mixing angle in matter 9 m . Results for both the infinite exponential profile 
and the BP2000 profile are shown. Fig. 0(B) shows the graph of A9 m = 
arcsin \ip2{9 m )\ — m as a function of 9 m . The main contribution to A9 m 
comes from the region 1.4 < 9 m < it/ 2, where the profile is close to the 
exponential. 

Fig. H shows the same evolution in the physical x-space. One can see that 
most of the contribution to A8 m comes from the part of the profile between 
0.7 'R Q and R & . This confirms that it would be incorrect to try to estimate 
A9 m by computing the slope around the point A = A, since for large 9 the 
entire region x > 0.7R Q is important. 

4 All known solutions can be derived from 
the one for the tanh profile. 

in Clon+mn 10 II TTri + V, 

m 



The formulation of the evolution equations introduced in Section |2.1| , with 9, 
as an independent variable, makes it possible to uncover a simple relationship 
between the known analytical solutions. Such solutions have been obtained 
for a very limited set of profiles. In addition to the exponential density 
distribution, explicit formulas in the literature have only been given for (i) a 
linear density distribution A(x) = —Cqx []22], |23|[] , 



(P c ) (lin) = exp(-7TA 2 (C )- 1 sin 2 29), (4.1) 

(ii) a distribution A(x) = B /x |35|] , 

exp(47rff cos 2 fl) - 1 
exp(47riJ ) — 1 

and (iii) the hyperbolic tangent distribution |3^] A(x) = A [l — tanh(x//)]/2, 
= cosh(7rM ) - cosh(7r/(A 00 - A)) 

{ cj(tanh) - cosh(7r/(Aoo + A)) _ CO sh(7r/(A 00 - A)) ' 1 J 



*Notice that Eq. (4J) is often written in a form containing a factor of 1/ cos 29. This is 
done to express Cq 1 in terms of the logarithmic derivative of the density at the resonance. 
In that form, it may superficially appear that P c has a singularity at 9 — tt/4. Of course, 
once the derivative is computed, the factors of cos 29 cancel out. As we have discussed 
in this paper, the conventional definition of the resonance has little physical meaning for 
large angles, so there is no good reason for introducing the factor of 1/ cos 29. 
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In the last equation A ro is the value of A m deep inside matter, A 



lim^oo A m = Al - 2A Q A cos 20 + A 2 . 

What do these distributions have in common that makes them exactly 
solvable? One important feature that unites them is that the corresponding 
differential equations can all be put in the hypergeometric form |J7|] . We next 
show more directly that Eqs. ( pT28|) , Q), (|4~2D , and (fO|) are all related to 



each other and that the first three can be easily obtained from the last one. 

To make this relationship clear, it is useful to show the analogs of Eq. 
( 2.15|) for the other three distributions in question. A straightforward calcu- 



lation yields 

/A m \ 2A 2 sin 2 2# 



(4.4) 



for the linear distribution, 



2B n sin 2 29 



\0m)(x/ x) sin 2 °m sin 2 (29 m - 29) 
for the 1/x distribution, and 

/ A m \ A l sin 29 sin 20 ^ 



(4.5) 



(tanh) 



sin 20 m sin(2# m - 20) sin(2# m - 20 o 



(4.6) 



for the hyperbolic tangent distribution (0^ is the value of m deep inside 
matter, so that sin 20^ = Asin20/A oo ). 

By inspecting Eqs. Q2.15| ), Q4.4|) , fl4.5|) , and (|4.6|) one can clearly see the 



common origin of all four solutions. For the hyperbolic tangent distribution 
(A m /^ m )( tanh ) has three simple poles in m on (0, vr/2], thus representing the 
most general case of all four. The exponential case is obtained when two of 
the poles, #oo and tt/2, merge in one double pole, and the 1/x case is obtained 
when the poles and Qoo merge. Finally, if all three poles merge in one triple 
pole at (tt/2), one obtains the linear case. Thus, given the result for the 
hyperbolic tangent distribution, it should be possible to recover the answers 
for the other three distributions by simply taking appropriate limits. 

We first show how the exponential result can be obtained from the one 
for the hyperbolic tangent. By comparing Eqs. ( |2.28| ) and (|4.3|) we see that 



25 



we need a large A limit of ( [4.3| ), since in this caseQ 9^ — > 7r/2, and also a 
substitution I — > 2r . 



ej (tanh) 



exp(vr/A )/2 - exp(7r/(A oc - A))/2 
exp(7r/(A 00 + A))/2 - exp(7r/(A 00 - A))/2 
exp(7r/(^o - (Aqq - A))) - 1 

exp(7r/(2A)) - 1 
exp(nl(A - A)(l - (A/A ) cos 20) + A)) - 1 
" exp(7r/(2A)) - 1 

exp(vr/A(l + cos20)) - 1 
exp(7d(2A)) - 1 
_ exp(47rr A cos 2 9) - 1 
exp(47rr A) — 1 

Above we used the fact that for large A A m = \J A$ — 2A A cos 29 + A 2 — ► 
Aq(1 — (A/Aq) cos 20). The physical interpretation of this result is the follow- 
ing. For a fixed A and A — > oo the part of the neutrino trajectory where adi- 
abaticity is maximally violated occurs at large x, where Ao[l+tanh(a;/Z)]/2 — > 
A exp(— 2x/l) = A exp(— x/r ). Notice, that A dropped out and the 
derivation is valid for all A and 9. 

To obtain the expression for P c for the linear profile from that for the 
exponential, in Eq. ( [2.15| ) we take the limit r —>■ oo, 9 — > 0, such that the 
product 2Arosin 2 26 l approaches a constant value, 2(A 2 /Co) sin 2 29. In this 
limit exp(47rroA cos 2 9) 3> 1 and sin 2 29 — > 4 sin 2 9, so that 



(-Pc)exp = a\ 1 * exp(-47rr Asin 



exp(47rr A cos 2 

exp(47rr A) — 1 
exp(-7r(A 2 /C ) sin 2 29) = (P c ) lin . (4.8) 



At last, we will show how the result for the 1/x distribution follows from 
that for the hyperbolic tangent distribution. The logic is similar to what was 
done before: Eq. ( |4.5| ) can be obtained from Eq. ( |4.6| ) by sending 9^ — ► 9 m 
and relabeling A l —+ 2A ; the result for the 1/x profile should then be read 
off from the solution of the differential equation with the profile (4.6|). A 



complication arises because in this case one needs to know the solution of the 



^Notice that the product Aosin2^oo in the limit of large Aq approaches a constant 
value A sin 29. 
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differential equation between ir/2 and 9^, while the known result, Eq. 
describes the solution between 6*00 and 9. The range [6*00, n/2] (and also [0, 9}) 
correspond to a different matter distribution, A(x) = Aq[1 + 1/ tanh(x//)]/2 
(see Appendix [A]). 

This difficulty can be easily resolved and the result in question can be 
obtained from Eq. (|4.3| ) by appropriate substitutions. The key observation is 
that in the case of the hyperbolic tangent the differential equation is uniquely 
specified by the relative positions of the three poles of Eq. ( [Of ) and the factor 
in the numerator. By shifting the poles such that 9 — > 0, we can reduce the 
problem of finding the solution on [#00,71-/2] to the solved case [9',9' OQ \ and 
use Eq. Q. 

The details of this procedure can be found in Appendix [A| After a 
straightforward calculation one finds that 

K = 



v sin(26»oo-26») 



^sh (Mo (-(ff^y + l)) - cosh (nlA (-^_ 

cosh (tcIA) — cosh (^(A^ — A )) 
cosh (7r/(Aoo + Aq)) — cosh (jrl(Aoo — A )) 



(4.9) 



From the first form of the expression it is easy to see that the formula satisfies 
the necessary nonadiabatic limit: lim^o P' c = cos 2 9. From the second form 
it is easy to take the limit which reproduces the formula for the 1/x profile. 
We can achieve 9^ — > 9 by making A large while keeping A fixed. Once 
again, in this limit Aoo — » A — A cos 29 and so 



P' 

c 



exp(7r/A) — exp(7r/(A 00 — Aq)) 
exp(7r/(A 00 + A Q )) - exp(7r/(A 00 - A )) 

exp(7T/(A - Aqq + Aq)) - 1 

exp(27r/A ) - f 
exp(27r/A cos 2 6') - 1 
exp(27r/A ) - 1 



(4.10) 



Upon relabeling A l — ► 2B we recover Eq. ( |4.2|) . 

The last result deserves a few comments. It is quite remarkable that the 
dependence of the jumping probability on the mass-squared splitting and 
energy completely dropped out at the end. The 1/x profile thus represents a 
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rather unique case when the adiabaticity is entirely determined by the density 
profile and the mixing angle. This could have been, of course, anticipated 
already on the basis of the dimensional analysis. Indeed, the parameter Bq is 
dimensionless and thus, together with 9, completely determines the jumping 
probability. 

It is instructive to see how the cancellation happens for small 9 and 
Bq 3> 1. Eq. (|4.2j ) then becomes P c = exp(—A-n-B 9 2 ). But we can also 
compute this differently: since for small angles the resonance is narrow, we 
can use a linear formula P c = exp(— nA 2 \dA(x)/dx\~ e l sin 2 29). Since for 
A(x) = B Q /x we have dA(x)/dx = —A(x) 2 /B and at the point of resonance 
A(x) = A, the quantity A cancels out of the final result: 

P c = exp(-nA 2 A9 2 B /A 2 ) = exp(-4ir9 2 B ) (4.11) 

The slope at the resonance point changes exactly in such a way as to com- 
pensate for the change in A 2 . The solar neutrino fits would be qualitatively 
different, were the solar density profile close to 1/x instead of the exponential. 

5 Conclusions 

In this paper we have discussed several aspects of the neutrino evolution in 
matter, emphasizing the case of the large values of the mixing angle, includ- 
ing values greater than 7r/4. We have pointed out how some of the results 
originally derived for small mixing angles can be modified to be applicable 
to all values of 9. Such results include the adiabaticity condition and the 
role played by the resonance in determining where the nonadiabatic jumping 
between the states takes place. We have formulated analytical criteria for 
the case of the exponential matter distribution and commented on how these 
criteria apply to case of the realistic solar profile. Although the focus of our 
analysis was on solar neutrinos, the results are useful for understanding the 
physics of neutrino oscillations in matter in general. 

We have presented the results of accurate numerical calculations, show- 
ing how the jumping probability P c interpolates between the QVO and the 
standard MSW regimes in the case of the realistic solar profile. An empirical 
prescription on how to estimate P c anywhere in this range with only elemen- 
tary functions was given. The matter effects in the quasivacuum regime were 
discussed. 
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Finally, we have shown that the known analytical solutions for the linear, 
exponential, and 1/x density distributions can be easily obtained from the 
result for the hyperbolic tangent. It was especially easy to see this using 9 m 
as an independent variable. In the process of the proof we also obtained an 
answer for a fifth distribution, N e oc (coth(x//) ± 1). 
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A Derivation of the expression for P c for the 
density distribution A(x) oc (coth(x//) ± 1] 



As mentioned in Section [|, one of the four matter distributions for which 
exact expressions for P c have been obtained is the profile A(x) = Aq{1 + 
tanh(x//)]/2. The answer is given by Eq. ( |4.3|) and represents the result of 
solving the differential equation ( |2.11|) on the interval [0,#oo], with /S. m /9 n 



given by Eq. ( f4.6|) . In this Appendix we show how this result can be used to 
obtain the solution to the same differential equation on the interval [0^, 7r/2]. 

Clearly for any x in the distribution A(x) = A [l + tanh(x//)]/2 the angle 
9 m is constrained between 9 and #00. The range of 9 m [6 l 00 ,7r/2] corresponds 
to a different density profile, which everywhere satisfies A(x) > Aq. This 
new profile can be found by simple analytical continuation. In the original 
profile the value of the density A(x) = C occurs at 

X=l 2 l ° g A^C- (A ' 1} 
When C > A the argument of the logarithm becomes negative. Using the 
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analytical continuation of the logarithm, it can be interpreted as 

X= 2K 0g \^C\ +m ) =X+ T- (A ' 2) 
Substituting this in the equation for A(x) we find 

A(x) = Aq[1 + tanh(x/Z + m/2)]/2 = A [l + 1/ tanh(x//)]/2. (A.3) 



To obtain the expression for P c for the distribution ( |A.3| ), as a first step 



in Eq. (|4.3|) we express the quantities Aoo and A m in terms of the angles 9 
and 6L 



'mi 



cosh(nlA Q ) - cosh (nlA ^%\ 
(tanh) «»h (vrMo-^^) - cosh (nlA ^^" * " ' ^ ' 



Next we shift 9 m — > 9 m — 9 in Eq. ( |4.6|) such that the pole at 9 moves to 0. 
This reduces the problem to solving the differential equation (|4.6|) between 
9'qq = Ti / 2 — 9 and 9' = 9^ — 9. To complete the change to the primed 
variables, in the numerator of Eq. ( |4.6| ) we substitute A sin 29 sin 29^ by 
A' sin 29' sin 29'^, where A' Q = A Q sin 29^/ sinffl^ - 29). We then find 



P' = 

c 

cosh(7r/A[ ] ) - cosh {^lA' ^^-le^) 

cnsh (tt1A> sin2e ^+ sin2e, N | _ nosh fvrM' sin2 ^° 
cosn [TTL/i sin(2dloB _ 2dl) ) cosn [ni/i sin(2e . 



-sin 26*' 



-20') 

^sh (^o^^) - cosh (vr/Ao ( ; - * 



v sin(26»oo -29) 



^sh (tt/Ao (-(ig^y + l)) - cosh (vr/A (^^y - 1 



(A.5) 



This is the answer that is needed in Sect. [|. 

One may be interested in the expression for P c for a distribution A(x) = 
A [l/ tanh(x/Z) — 1]/2, which has the property that it vanishes in the limit of 
large x. This can be found by redefining the vacuum values of the neutrino 
oscillation parameters to be A^ and 9^. Eq. ( |A.5| ) can then be rewritten as 



, _ cosh (ttZA) - cosh (^/(A^ - A )) 

c ~ cosh (vr/( Aoo + A Q )) - cosh (ttZ ( A^ - Aq)) ' 



where A = J A% + 2/loAoo cos 29^ + A^. 
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